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Abstract 

We present a short review of the existing evidence in favor of neutrino mass and 
neutrino oscillations which come from different kinds of experiments. We focus our 
attention in particular on solar neutrinos, presenting a review of some recent anal¬ 
ysis of all available neutrino oscillation evidence in Solar experiments including the 
recent SNO CC and NC data. We present in detail the power of the reactor experi¬ 
ment KamLAND for discriminating existing solutions to the SNP and giving accurate 
information on neutrino masses and mixing angles. 
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1 Introduction 


Neutrino has a long and fascinating history that started more than seventy years ago, 
when Pauli postulated its existence [1] to explain the data about the high energy part of 
the (3 decay spectrum. Since then the main issue was to discriminate wether neutrino is a 
massive or massless particle and eventually to find the value of its mass. 

In this work we briefly review the main steps that have been done towards the solution 
of the neutrino mass puzzle and the different kind of experiments that have been studying 
it. Then we remind and critically discuss the essential data that have been obtained in 
this last exciting period and the ones that should be presented very soon, focusing in 
particular our attention on solar and reactor neutrino physics. We also present the salient 
aspects of the phenomenological analysis of these data that have been developed by our 
group, with the aim to understand at which level of accuracy we can determine the values 
of the mixing parameters and how we can use the results of the forthcoming experiments 
to improve this accuracy. 

The first studies of neutrino mass, performed in the 40th [2], were based on the so 
called Fermi-Perrin method of observation of the (3 spectrum near the end point and they 
gave the limit m v < 500 eV. This limit has been significantly lowered in the following 
decades and by now we know, for instance, that m{y e ) < 2.2eV. About ten years later 
than the first measurement of neutrino mass, Goldhaber in ’58 measured its helicity [3] 
and proved that the neutrinos pruduced in (3 decays behave like left-handed particles. 

An essential idea for the understanding of neutrino properties was the oscillation 
hypothesis introduced by Pontecorvo in ’57 [10]. He conjectured that the neutrino states 
produced in weak interactions could be superpositions of two different states of Majorana 
neutrinos, having different values of the mass; therefore there could be an oscillation from 
one neutrino state into another, in a way very similar to what happens in a neutral kaon 
system. This revolutionary idea (implying the existence of at least one mass eigenvalue 
different from zero for neutrinos) has been further developed in a series of works [11] and 
it has been confirmed by many experimental evidences in the following years, as we will 
see in the rest of this paper. Therefore, by now, we are almost sure that neutrinos are 
massive and oscillating particles. 

The experimental hints in favor of a non zero neutrino mass cannot be accomodated in 
the usual “minimal” version of the Standard Model of electroweak interactions, in which 
neutrino is simply a left-handed particle and it is impossible to build a renormalizable 
neutrino mass term. On the other hand, the existence of such a mass term is predicted 
in most of the theories beyond the Standard Model, like, for instance, supersymmetric 
and unified theories. Even if the main issue whether neutrino is massive or not seems 
to be solved, there are still essential theoretical questions to answer. We don’t know 
yet, for instance, whether neutrino is a Dirac or Majorana particle, neither we have an 
unique explanation of the reason for which it is so lighter than all the other known massive 
particles. The mechanism usually invoked to explain neutrino lightness is the so called see¬ 
saw mechanism [4], but there are also other possible alternative explanations. We don’t 
know yet the absolute value of neutrino masses (experiments up to now can only put upper 
limits and give information on the mass differences) and which is the exact structure of 
neutrino mass spectrum. We are not able, for instance, to decide whether we are in 
presence of “normal” or “inverted hierarchy” or if the masses are “quasi degenerate” [5]. 

For all these reasons, we can say that in the following years neutrino physics will 
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continue to be an ideal playground for testing different theories beyond the Standard 
Model, like, for instance, the supersymmetric and unified theories [6] or the ones predicting 
the existence of large extradimensions [12]. This last class of models, which become very 
popular in the last years, have important consequences also for neutrino phenomenology. 
For instance, they predict the existence of an infinite tower of sterile neutrinos and they can 
give rise to mechanism explaining the lightness of neutrino mass, similarly to the see-saw, 
but without the need of having a very high mass scale. However the recent experimental 
results, mainly from cosmology, put very stringent constraints on the large extradimension 
theories and ruled out the simpler versions of these models. 1 . 

The experiments aiming to measure neutrino mass or at least to find evidences that 
it is different from zero can be classified in three big categories: 

• the direct kinematical searches; 

• the searches of neutrinoless double (3 decay; 

• the experiments looking for neutrino oscillations as a proof that neutrinos are mas¬ 
sive. 

The direct kinematical searches are performed by looking at the high energy part of 
the 3 H /3-decay spectrum. The present limits on the values of v T and masses are [8, 9]: 

m{y T ) < 18.2 MeV 

m{v lx ) < 190 KeV. (1) 

The best limits for the mass of the electron neutrino, instead, have been obtained from the 
Mainz and the Troitsk [23, 24] experiments that have found m(y e ) < 2.2 eV. In future 
many experiments will try to lower this limit. In particular there is a great expectation 
for KATRIN (the Karlsruhe Tritium Neutrino experiment) [25], that should start the data 
taking in 2007 and improve the sensitivity down to 0.35 eV. 

The search for neutrinoless double j3 decays is very important because the observation 
of these decays (violating by two units the lepton number conservation) would imply, 
under the assumption of CPT invariance, the Majorana nature of neutrino [18,19]. The 
most stringent limits available at the moment for the effective Majorana mass come from 
the Heidelberg-Moskow collaboration [21] (ro„) < 0.35 eV and from IGEX (International 
Germanium Experiment) [22] {m u ) < 0.33 — 1.35eV. 

The effective Majorana mass is given by: 

|(m)| = |mi|t/ e i| 2 + m 2 |t/ e2 |V Q21 +m 3 |t/e 3 |V Q31 |, (2) 

where m, indicate the values of the mass eigenstates, the U e i are elements of the Pontecorvo- 
Maki-Nakagawa-Sakata mixing matrix and a 2 i,a 3 i are the two additional CP violating 
Majorana phases. The masses and the elements of the mixing matrix entering equation (2) 
can be expressed in terms of the mass differences and mixing angles that can be recov¬ 
ered from the experiments on solar, atmospheric and reactor neutrinos. Thus, given these 
experimental inputs, the value of |(m)| depends on the ones of m\ and of the Majorana 
phases and on the general structure of neutrino mass spectrum (i.e. normal or inverted 

1 About the recent status of the models predicting large extradimension and the relative experimental 
constraints see, for instance, the talk of J. Liu at this conference 
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hierarchy or quasi degenerate case). In particular one can get interesting lower bounds in 
the cases of inverted hierarchy and quasi degenerate mass spectrum and a stringent upper 
bound [45] (smaller than lCU 2 eV) in the case of normal hierarchy. The NEM03 [42] ex¬ 
periment and the criogenic detector CUORICINO [20] are expected to reach sensitivity of 
the order of 10 _1 eV and even an order of magnitude better is planned for the next genera¬ 
tion experiments [20,44] CUORE, GENIUS, EXO, MAJORANA and MOON. Therefore, 
there is the hope that the results of these future experiments, complemented with the 
information on the absolute value of the masses that could come from KATRIN and from 
cosmology, could be used to recover hints of CP violation in the leptonic sector or to 
discriminate between the different possible structures of the mass spectrum. For a more 
detailed analysis of the present situation and the future perspectives of neutrinoless double 
/3 decays we refer the interested reader to the work [46] and the references contained in it. 

A first group of experiments looking for oscillation signals use neutrino fluxes pro¬ 
duced at accelerators and nuclear reactors. They are usually distincted in long- and short 
baseline, according to the distance existing between the neutrino production point and 
the detector. Many short baseline accelerator experiments, like NOMAD [26] and CHO¬ 
RUS [27] at Cern, didn’t find any signal of oscillation and, consequently, they gave im¬ 
portant constraints on the possible values of the mixing parameters. Between the reactor 
experiments, it is worthwhile to remember the results of CHOOZ [28] and Palo Verde [30]. 
At CHOOZ a beam of reactor v e was sent to a detector located about 1 Km far away and 
it was detected through the reaction v e +p e + + n. No evidence of oscillation was found 
and these results can be used to exclude a significant part of the mixing parameters plane. 
In particular they tell us that Am 2 must be smaller than 10 _3 eV 2 , unless the values of 
the mixing angle are very small. The CHOOZ data can be used also to fix an upper limit 
for the quantity sin 2 (6), where sin 2 6 coincides with | E/e 3 1 2 in the case of normal hierarchy 
and sin 2 6 = |U e i| 2 for neutrinos with inverted hierarchy mass spectrum. 

The opposite situation took place in the case of LSND [31], a short baseline accelerator 
experiment performed with a neutrino beam produced at the Los Alamos meson physics 
facility (LAMPF). This experiment found evidences for two kind of oscillation signals: 
that is an excess of v e in the beam of n lt produced by the decay at rest of the /j, + (obtained 
as secondary products of the proton accelerator beam), and a signal of oscillations into 
u e of the produced by the 7r + decay in flight. The LSND result, if confirmed, would 
be a clear indication of oscillation with very high values of the mass difference, up to 
A777 2 — leV 2 . To reconciliate this result with the ones coming from solar and atmospheric 
neutrinos (that we are going to present), one would have to postulate the existence of 
at least one sterile neutrino in addition to the usual three active ones. However, up to 
now there have been no independent confirmations of the LSND results. The KARMEN 
experiment [33], performed at the Rutherford Laboratories, explored a significant part of 
the mixing parameter space proposed by LSND and it didn’t find any signal of oscillation 
2 . A new experiment MiniBoone [35,36] just started running and it shuold soon produce 
data. It has been projected in such a way to test definitely the validity of LSND results. 

A new generation of very long baseline experiments became available in the last years. 
The forerunner of them is K2K [37, 38], that uses a neutrino beam produced at the Japan 
kaon facility KEK and detected at the Kamioka site. Up to now K2K has detected 56 
events instead of the expected value in absence of oscillations, that is 80lg' 3 . This is 

2 About the compatibility of LSND and Karmen results see also [34] 
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a confirmation of neutrino oscillations (the null oscillation probability is less than 1 %). 
Moreover the best fit point [38] values for the mass difference and the mixing angle (Am 2 = 
2.8 x 10~ 3 eV 2 and sin 2 26 = 1) are in good agreement with the results of atmospheric 
neutrino experiments. Two similar projects have been already approved and will become 
available in the near future: one of them is a neutrino beam from CERN to the Gran 
Sasso Labs [13] and the other one is in the States [14] (from FNAL to Soudan). The 
long baseline accelerator experiments will probably give an important confirmation of the 
oscillation evidences coming up to now from the study of solar and atmospheric neutrinos. 
They are also expected to find in an unambiguous way indications of oscillation from 
appearance signals. In the long baseline experiment one has in addition the opportunity 
of choosing the specific characteristic of the beam; hence they can be used to perform 
precision measurements. For instance they should be useful to study the value of the 
mixing angle @ 13 , relevant for eventual CP violation. The present limit on the measurement 
of this angle coming from CHOOZ ( 6*13 < 9 degrees), could be lowered to the level of about 
5 degrees at ICARUS, one of the two experiments that will use the CERN-Gran Sasso 
beam. 

Important result should very soon come also from the long baseline reactor experi¬ 
ment kamLAND [48], that could in principle give a definite solution to the solar neutrino 
problem, as we are going to see in the following. 

The two main categories of experiments looking for oscillation signals are the ones 
that study the atmospheric and the solar neutrinos. 

The flux and the angular distribution of electronic and muonic neutrinos produced in 
the atmosphere, in the processes generated by cosmic rays, is known with quite a good 
accuracy [15], taking into account the number and the properties of cosmic rays and the 
important geomagnetic effects. This input information are essential, because most of the 
atmospheric neutrino experiments measure the value of the double ratio R = . 

The numerator and denominator are, respectively, the experimental and the Monte Carlo 
computed values of the ratio between the events generated by muonic neutrinos (and 
antineutrinos) and the ones generated by electronic neutrinos (antineutrinos). 

At present there are different experiments running on atmospheric neutrinos and 
most of them are essentially water Cherenkov detectors (like Kamiokande [49-52], Super- 
Kamiokande [53-56], IMB [57]) or iron plate calorimeters (like Soudan II [17] and in the 
past years Frejus [59] and Nusex [61]). Clear evidences of oscillation have been found at 
Kamiokande, Super-Kamiokande (SK), IMB and Soudan II and also at the MACRO [16] 
experiment at Gran Sasso. 

Here we recall the results of SK, which are based on a very high statystic and are, 
in any case, in good agreement with the other atmospheric experiments which found os¬ 
cillation signals. The SK results are R = 0.638 ± 0.016 ± 0.050 for the Sub-GeV events 
and R = 0.658lg'Q 3 g ± 0.078 for the Multi-GeV events. Another interesting observable is 
the up-down asymmetry between the up going events, in which the neutrino crossed the 
Earth before interacting in the detector, and the down going ones: A e ^ = The 

experimental value of this quantity is consistent with zero for the electronic neutrinos, 
while for the muonic ones the day-night asymmetry for high values of the momenta is a 
decreasing negative value. This result is a clear indication of a reduction of the flux of 
muonic neutrinos and antineurinos that arrive to the detector after crossing the Earth. 
The most natural explanation of this phenomenon is the possibility that the muonic neu¬ 
trinos oscillate into other flavors and the oscillation probability is greatly enhanced by the 
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interaction with matter. 

The last group of experiments is that of the experiments observing the neutrinos 
coming from the Sun. We will discuss them in detail in the rest of the paper. 

1.1 History of the solar neutrino problem 

The first experiment on solar neutrinos, Homestake [65], started at the end of the ‘60s 
using the inverse (3 decay on chlorine 37 Cl + v e — A 7 Ar + e~. The threshold energy 
was E t f ir ~ 0.81MeV, hence it was sensitive to the pep, ' Be , 8 H and hep components of 
the solar neutrino flux. The results were really surprising, because Homestake found a 
deficit of the solar neutrino flux of more than 60% that predicted by the Solar Standard 
Model (SSM). The updated value of the ratio R, between the experimental results and 
the SSM prediction [66], for the chlorine experiment is R = 0.34 ±0.03. This result raised 
fundamental questions: what happens to solar v on their way to earth? Eventually, could 
the SSM be wrong? 

The Homestake indication was confirmed by similar experiments, SAGE [70] in Russia 
and GALLEX [71] and later on GNO [72] at the INFN Gran Sasso Labs, which used gallium 
instead of chlorine. The energy threshold is lower in the gallium experiments (E t h r ~ 
233 keV) making them also sensitive to pp neutrinos, which are the main component of 
the solar neutrino flux. The updated gallium results are [70-72]: 

R = 0.60 ±0.05 (SAGE) 

R = 0.58 ±0.05 (GALLEX - GNO) . (3) 

This confirmation of Homestake results gave a strong support to the neutrino oscillation 
hypothesis and caused an increase of the interest for this problem. It could be a signal of 
new physics. 

An essential improvement in the knowledge of solar neutrinos came with the advent of 
the water Cherenkov experiments, Kamiokande [73] and Super-Kamiokande (SK) [74, 75], 
that looked at the elastic scattering v e ± e~ —> v e ± e~ and confirmed the existence of the 
“solar neutrino problem” with a very high statistic. In this experiments it was possible 
to know the direction of the incoming neutrino (by looking at the outgoing direction of 
the recoil electron) and also to study the energy and angular spectrum and the day night 
asymmetries. The energy threshold for these experiments was quite high (5 MeVfor SK) 
and therefore they were sensitive only to the high energy component of the neutrino flux, 
that is 8 1? and hep neutrinos. Their results confirmed the existence of a deficit in the 
electron neutrinos reaching the detector. The SK result for the energy spectrum and 
the small values of the day-night asymmetries were also very important to put strong 
constraints on the possible values of the mixing parameters. 

After the publication of SK data it was clear that there was a deficit of solar electron 
neutrinos reaching the Earth, with respect to the flux predicted by SSM. The oscillation 
hypothesis was considered the most plausible explanation of this phenomenon, but there 
were still different regions allowed by the experiments in the mixing parameter plane, as 
we will see in detail. 
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1.2 The post SNO situation 

The real breakthrough was due to the SNO experiment that published its first data in 
2001 [76] ( see also [82]). SNO is a deuterium Cherenkov detector designed to look simul¬ 
taneously at three different reactions: 

e~ + p + p (Charged Current)-, 
e~+n + p (Neutral Current)-, (4) 

v x + e~ (Elastic Scattering). 

The first reaction (CC) receives contribution only from the electron neutrino, while the 
others (NC and ES) are sensitive to all neutrino flavors. This experiment gives the first 
direct model independent measurement of the total solar neutrino flux reaching the Earth 
(through the NC observation) and at the same time, comparing this flux with the one of 
v e recovered from CC, it offers a strong evidence of the oscillation of v e into other active 
neutrinos. 

During its first phase of working [76] SNO observed the charged current and elastic 
scattering events, with an energy threshold for electron detection of 6.75 MeV. The v e flux 
measured from CC, after 241 days of running, was: = 1.75±0.07(siai.)lo'ii(sysi.) x 

10 6 cm~ 2 s~ l . The ratio between this value and the SSM prediction was R = 0.35 ± 0.03. 
In the SNO experiment the neutrino flux can be recovered also from the elastic scattering, 
using the relation 3 : 

^ S = <S>l s + 0.154 Y, ( 5 ) 

j=/i,T 

The value of the total neutrino flux recovered from the elastic scattering at SNO and also, 
with a better statistics, at SK doesn’t agree with the v e flux obtained from SNO CC. The 
comparison of the two results gives: 

Y 4>f; S = 3.69 ± 1.13 x 10 6 crrUV 1 . (6) 

i=H,r 

This result was the first evidence (at 3cr level) of the presence of muonic and tauonic 
neutrinos in a electronic neutrino beam reaching the Earth from the Sun. Therefore it 
was, up to the present SNO data on NC, the most robust evidence of v e oscillation into 
other active neutrinos. It is also remarkable that the sum of the v e and u thT fluxes give a 
value in good agreement with the SSM prediction. Consequently the results of SNO phase 
I also strongly disfavored the hypothesis of pure oscillation into sterile neutrinos. 

Observation of neutral current neutrino interactions on deuterium in the SNO ex¬ 
periment has been recently presented [77,78]. Using the NC, ES and CC reactions and 
assuming the 8 B neutrino shape predicted by the SSM, the electron and active non-electron 
neutrino component of the solar flux at high energies ( > 5 MeV) are obtained. The non¬ 
electron component is found to be ~ 5 a greater than zero, the standard prediction, thus 
providing the strongest evidence so far for flavour oscillation in the neutral lepton sector: 
the agreement of the total flux, provided by the NC measurement with the expectations 
implies as a by-product the confirmation of the validity of the SSM [67-69]. 

3 the contribution of and v T to the elastic scattering cross section is only through neutral current, 
hence it is about 1/6 of the contribution of ry, that can interact also through charged current 


V e + d 

v x + d 

Vx + e' 


7 



The results presented recently by SNO on Solar neutrinos [84] confirm and are consis¬ 
tent with previous evidence from SK and the rest solar neutrino experiments [85, 86, 96, 97]. 
The CC, ES and NC global and day and night fluxes presented in Refs. [77, 78]. are derived 
under the assumption that the S B spectral shape is not distorted from the SSM prediction. 
With this assumption the SNO collaboration checks the hypothesis of non-oscillation, or 


zero 4>[j,+t flux. 

Different quantities can be defined [80] in order to make the evidence for disappearance 
and appearance of the neutrino flavours explicit. Letting alone the SNO data, from the 
three fluxes measured by SNO is possible to define two useful ratios, deviations of these 
ratios with respect to their standard value are powerful tests for occurrence of new physics. 
Here we compute the values for (frcc/^ES and (^Cc/^NC being specially careful with the 
treatment of the correlations on the uncertitudes, the inclusion or not of these correlations 
can affect significantly the results for these ratios (see table II in Ref. [77]. for a complete 
list of sistematical errors). For the first ratio, from the value from SNO rates [77, 78]. we 
obtain 


<ftcc 

4>es 


n 73+0.10 


a value which is ~ 2.7 a away from the no-oscillation expectation value of one. The ratio of 
CC and NC fluxes gives the fraction of electron neutrinos remaining in the solar neutrino 
beam, our value is: 


4>CC 

0NC 


n 34+0-05 
u -'+ i -0.04> 


this value is nominally many standard deviations (~ 13a) away from the standard model 
case [41]. 

Finally, if in addition to SNO data we consider the flux predicted by the solar standard 
mode one can define, following Ref. [98], the quantity sin 2 a, the fraction of ” oscillation 
neutrinos which oscillated into active ones”, again using the SNO data and fully applying 
systematic correlations (see table 2 in Ref. [77, 78]), we find the following result: 


. 2 4>nc - 4>CC 

sin a = --;- 

4>SSM ~ <PCC 


0.92 


+0.39 

- 0 . 20 - 


The SSM flux is taken as the S B flux predicted in Ref. [69]. Note that, although consistent 
with it, this result differs significantly from the number obtained in Ref. [98], this is due 
to the introduction of systematic correlations in our calculation. The central value is 
clearly below one (only-active oscillations): the fraction of sterile neutrinos is cos 2 a < 0.28 
(la). Although electron neutrinos are still allowed to oscillate into sterile neutrinos the 
hypothesis of transitions to only sterile neutrinos is rejected at nearly 5a, this significance 
would be reduced if we consider applying a 1-sided analysis to avoid non-physical values. 


2 Global analysis of the solar neutrino data 

Given all the experimental data that we have just reported, one can say that there is really 
strong evidence that neutrinos are massive and oscillating particles. Nevertheless, details 
of the mass pattern still have a long way before being clarified. With this aim in mind, 
we developed a global analysis of all the available data on solar neutrinos, also including 
the CHOOZ constraints. Our first purpose was that of determining the regions in the 



mixing parameter plane that are still compatible with the experiments. In addition to 
this, we wanted to understand how the forthcoming experiments (in particular Borexino 
and KamLAND) could improve our knowledge of neutrino mass properties. 

We assumed neutrino oscillation as a working hypothesis and considered bidimensional 
models. For details of our analysis we refer the interested reader to [79, 80]. Here we just 
report the most salient aspects of our strategy. The analysis is based on the numerical 
calculation of the expected event rate for every solar neutrino experiment as a function of 
the mixing parameters and on the comparison between these expected numbers and the 
experimental data. The statistical analysis is based on the y 2 method its output being 
contour plots in which one can see which values of the mixing angles and mass differences 
are still allowed at a given confidence level. 

Our calculation can essentially be split into two parts. The first one is the determi¬ 
nation of the neutrino transition amplitude, i.e. the probability for an electronic neutrino 
produced in the Sun to change its flavor before reaching the detector. The other ingredi¬ 
ent is the calculation of the detector response functions, that, for a given neutrino energy, 
depend on the experimental details of the specific detector (i.e. efficiency, resolution, etc.) 
and on the cross section for the reaction under examination. 

The transition amplitude calculation is separated in three parts, corresponding to the 
neutrino propagation inside the Sun, in the vacuum and in the Earth. For every value 
of the mixing parameters we compute fully numerically the amplitudes in the Sun and in 
the Earth, while the one corresponding to the vacuum evolution is computed analytically. 
The three amplitudes are patched together using the evolution operator formalism [81]. 

3 The present situation 

We included in our analysis [80] the total rates of the chlorine and gallium experiments, 
together with the different energy bins of SuperKamiokande and with the most recent 
results of SNO [77, 78]. The results are shown in Fig.l where we have generated acceptance 
contours in Am 2 and tan 2 $. In Table (1) we present the best fit parameters or local 
minima obtained from the minimisation of the y 2 function. Also shown are the values 
of Xmin P er degree of freedom (x 2 /n) and the goodness of fit (g.o.f.) or significance level 
of each point (definition of SL as in Ref. [9]). In order to obtain concrete values for the 
individual oscillation parameters and estimates for their uncertainties, it is preferable to 
study the marginalized parameter constraints. It is justified to convert \' 2 i n t° likelihood 
using the expression C = e~ x / 2 , this normalised marginal likelihood is plotted in Figs. (2) 
for each of the oscillation parameters Am 2 and tan 2 6. For tan 2 6 we observe that the 
likelihood function is concentrated in a region 0.2 < tan 2 0 < 1 with a clear maximum at 
tan 2 6 0.4. The situation for Am 2 is similar. Values for the parameters are extracted 

by fitting one- or two-sided Gaussian distributions to any of the peaks (fits not showed in 
the plots). In the case of the angle distribution the goodness of fit of the Gaussian fit is 
excellent (g.o.f > 99.9%) even at far tail distances thus justifying the consistency of the 
procedure. The goodness of Gaussian fit to the distribution in squared mass, although 
somewhat smaller, is still good. The values for the parameters appear in Table 1. They 
are fully consistent and very similar to the values obtained from simple x 2 minimisation. 

In summary, the direct measurement via the NC reaction on deuterium of 8 B neutrinos 
combined with the CC results have largely confirmed the neutrino oscillation hypothesis. 
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Method 


A) Minimum LMA Am 2 = 5.44 x 10 5 eV 2 

B) From Fit Am 2 = 4.5l 2 x 10 _5 eV' 


tan 2 9 = 0.40 X 2 m = 30.8 g.o.f.: 80% 

tan 2 9 — 0.40^o og 


Table 1: Mixing parameters: A) Best fit in the LMA region and B) from fit to marginal likelihood 
distributions. 


10“ 3 


10“ 4 

> 

-it 

% 

<1 
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10“ 6 

10“ 2 10 -1 1 10 
tan 2 6 

Figure 1: Allowed areas in the two neutrino parameter space. The point with error bars corre¬ 
sponds to best the results from fit to marginal likelihoods. The coloured areas are the allowed 
regions at 90, 95, 99 and 99.7% CL relative to the absolute minimum. The region above the upper 
thick line is excluded by the reactor experiments [29]. 



SK Sp(d)+Sp(n) 

A w =7.0%± 5.1% 
0rB)=5.O9± 0.64 


We have obtained the allowed area in parameter space and individual values for Am 2 and 
tan 2 9 with error estimation from the analysis of marginal likelihoods. 

4 Future scenarios: Borexino 

Given this situation, we studied which new information should come in future from the 
Borexino data [79,83,87]. Borexino [88,89] is a solar neutrino experiment, mainly sensi¬ 
tive to the "Be component of the neutrino flux, that should start running in very next 
years at the Gran Sasso Labs. In Figure 3 the usual contour plots obtained from all the 
experiments available up to now are superimposed to the contour lines corresponding to 
different hypothetical possible values of the total rate at Borexino. Note that, for what 
concerns the data coming from the other solar neutrino experiments, this figure refers to 
the situation how it was in 2001, before the publication of the last results from SK and 
from SNO neutral current. As one can see from the picture, Borexino should be able to 
clarify the situation in the case in which the solution very well is in the small mixing 
angle region. The situation would be, instead, more complicate in case of LMA or LOW 
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solutions. In these two regions, in fact, the ratio between the Borexino signal and the 
SSM prediction in absence of oscillations should be between 0.6 and 0.7 . The discrimina¬ 
tion power of Borexino increases a lot if we look also at the day-night asymmetry, as one 
can see from Figure 4. The LOW region is characterized by high values of the asymme¬ 
try, that can reach up to 20%, while in the LMA region the day-night asymmetry is much 
lower. Hence, by looking simultaneously at the total rate and at the day-night asymmetry, 
Borexino should be able to discriminate between the two solutions of the solar neutrino 
problem that are compatible with the experiments up to now, that is the LMA and the 
LOW solutions. 

5 Future Experiments: Kamland 

Another very important experiment, already running, that should significantly improve 
our knowledge of the mixing parameters relevant for solar neutrinos is KamLAND [48]. 
In this experiment a flux of low energy v e produced by different nuclear reactors is sent 
to a scintillator detector capable of detecting their interactions with protons. Although 
it not a solar neutrino experiment, KamLAND is sensitive to neutrino oscillations with 
mixing parameters in the LMA region, that seems to be the solution of the solar neutrino 
problem preferred by the present data. Therefore, we can hope that KamLAND will soon 
be able to determine the exact values of the mixing parameters with satisfactory accuracy. 
The main limitation of KamLAND is its reduced sensitivity to the extreme upper part of 
the LMA region, that could create problems in the determination of Am] 2 , as discussed 
in [90,91] and later on in [92]. For a detailed discussion about KamLAND potentiality 
and discrimination power we refer the interested reader to [83,93]. 

In order to study the potentiality of KamLAND for resolving the neutrino oscillation 
parameter space, we have developed two kind of analysis. In the first case (Analysis 
A below) we deal with the KamLAND expected global signal. We asumme that the 
experiment measure a certain global signal with given statistical and systematic error after 
some period of data taking (1 or 3 yrs) and perform a complete y 2 statistical analysis 
including in addition the up-date solar evidence. In the second case, Analysis B, we 
include the full KamLAND spectrum information. Instead of giving arbitrary values to 
the different bins, we assume a number of oscillation models characterized by their mixing 
parameters (Am 2 , 9). After including the solar evidence we perform the same y 2 analysis 
as before. 

In Fig.4 we graphically show the results of this analysis. They represent exclusion plots 
including KamLAND global rates, given a hypothetical experimental global signal ratio: 
respectively strong and medium suppression S/So = 0.3,0.6 and no oscillation evidence 
S/So = 1.0 for one and three years of KamLAND data taking. As can be seen in the 
figures, as the KamLAND experimental signal decreases, the LMA region is singled out. 
The periodic shape in Am 2 of the 90 % C.L. (red regions) which becomes apparent in the 
three-year plot is due to the periodicity of the response function: in order to distinguish 
among these different equally-likely solutions, one must analyze the energy spectrum, 
Obviously, only if KamLAND sees some oscillation signal (i.e. Si/So « 1.0 ) does the 
LMA region become the only solution. If we consider a hypothetical signal closer to 1.0 
than 0.3, we see that the LOW region survives, although it is less favored. 

In Figs. (5) we graphically show the results of this analysis for a selection of points and 
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tan 2 6 Am 2 /eV 2 


Figure 2: Marginalized likelihood distributions for each of the oscillation parameters Am 2 (right), 
tan 2 9 (left). The curves are in arbitrary units with normalization to the maximum height. Values 
for the peak position are obtained by fitting two-sided gaussian distrbutions (not showed in the 
plot). Dashed lines delimit ±1 <t error regions around the maximum. 



tan 2 8 


Figure 3: Possible values of the total signal at Borexino. The full lines correspond to a total rate 
equal respectively to 0.5, 0.6 and 0.7 with respect to the SSM prediction in absence of oscillation. 
These isosignal lines are superimposed to the contour plots (colored regions) corresponding to the 
regions allowed at different confidence levels from the other experiments. 


12 

















for three years of data taking restricting ourselves to the LMA region of the parameter 
space where, as we have noted before, the KamLAND spectrum information is specially 
significative. In each plot allowed regions corresponding to different starting points are 
superimposed, every region is distinguished with a label. The position of initial points is 
labeled with solid stars. 

The first case, study of the KL spectrum alone (x 2 pec kl) represented by the Fig. (5 
left). The allowed parameter space corresponding to each particular point is formed by 
a number of, highly degenerated, disconnected regions symmetric with respect the line 
tan# = 1. These regions can extend very far from the initial point specially in terms of 
Am 2 but also in some occasions in terms of tan 2 9. For example the point “A” located at 
(Am 2 = 5.7 x 10~ 4 , tan 2 0 = 0.38) gives rise to two sets of thin regions situated respectively 
at Am 2 ~ 10~ 3 ,10~ 4 and a third region situated at Am 2 ~ 10~ 5 which practically covers 
the full range tan 2 9 ~ 0.1 —10. A similar behavior is observed for point “B”. Of course this 
situation is not very favorable for the future phenomenologist trying to extract conclusions 
from the KamLAND data. A much comfortable situation is found for points nearer the 
center of the LMA region. Note how the regions corresponding to the points “D,E” and 
specially “F” only extend very gently around the initial location. 

The results of the full analysis are summarized in Fig.(5, right). The position of the 
minima of y 2 , marked in the plot with crosses, is practically identical to the position 
of the initial points except in some case where the difference is not significant anyway. 
The general effect of the inclusion of the solar evidence in the x 2 is the breaking of the 
symmetry in tan 2 9 as expected and the general reduction of the number of disconnected 
regions corresponding to each point. Note however that the point “A” still gives rise 
to a small allowed region situated nearly one order of magnitude smaller in Am 2 . The 
“B” region is shrinked near its initial location as happens to the rest of points. The 
conclusion to be drawed from these plots is that KamLAND together with the rest of 
solar experiments will be able to resolve the neutrino mixing parameters with a precision 
of 5 log Am 2 ~ ±0.1 practically everywhere. However, for values of Am 2 > 10 -4 the 
problem of the coexistence of multiple regions with similar statistical significance will still 
be present. 

5.1 KamLAND and solar antineutrinos 

We have also investigated [94] the possibility of detecting solar antineutrinos with the 
KamLAND experiment. These antineutrinos are predicted by spin-flavor solutions to the 
solar neutrino problem. As we saw before, the recent evidence from SNO shows that a) 
the neutrino oscillates, only around 34% of the initial solar neutrinos arrive at the Earth 
as electron neutrinos and b) the conversion is mainly into active neutrinos, however a non 
e,/U — r component is allowed: the fraction of oscillation into non-^ — r neutrinos is found 
to be cos 2 a = 0.08 !q' 2 q. This residual flux could include sterile neutrinos and/or the 
antineutrinos of the active flavors. 

The KamLAND experiment is potentially sensitive to antineutrinos coming from solar 
8 B neutrinos. Reactor and solar antineutrino signals are shown in Fig.(5.1). In case of 
negative results, we find that the results of the KamLAND experiment could put strict 
limits on the flux of solar antineutrinos <1>( S R) < 1.0 x 10 4 cm~ 2 s -1 , and their appearance 
probability (P < 0.2 — 0.15%), respectively after 1-3 years of operation. Assuming a 
concrete model for antineutrino production by spin-flavor precession in the convective 
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tan 2 9 


Figure 4: Exclusion plots including KamLAND global rates (Analysis A), given a hypothetical 
experimental global signal ratio: respectively S/Sq = 0.3. Statistical and assumed systematics 
(~ 5%) errors are included. Three year of KamLAND data taking. The colored areas are allowed 
at 90, 95, 99 and 99.7% CL relative to the absolute minimum. The region above the upper thick 
line is excluded by the reactor experiments [29]. 



Figure 5: Allowed areas in the two neutrino parameter space after 3 years of data taking in Kam¬ 
LAND (Analysis B). Allowed regions belonging to the same point are labeled with the correspond¬ 
ing letter, the position of the point itself is labeled with a solid star. The colored lines separate 
allowed regions at 90, 95, 99 and 99.7% CL relative to the absolute minimum. (Left) Results 
with the KamLAND spectrum alone. (Right) KamLAND spectrum plus solar (CL,GA,SK,SNO) 
evidence. Crosses are situated in the position of the y 2 minima. 


14 









^thr 

Ssun 

SRet 

Bckg. 

P (CL 95)% 

P (CL 99)% 

6 MeV 

616 

43 

70 

0.22 

0.23 

7 MeV 

500 

11 

65 

0.19 

0.20 

8 MeV 

366 

2 

60 

0.21 

0.23 


Table 2: Expected signals from solar antineutrinos after 3 years of data taking. Reactor 
antineutrino (no oscillation is assumed) and other background (correlated background) 
over the same period. The random coincidence background is supposed negligble above 
these energy thresholds. Upper limits on the antineutrino oscillation probability. 



Figure 6: The KamLAND expected positron spectra in absence of oscillations coming 
from reactor antineutrinos (normalized to three years of data taking). The “solar” positron 
spectrum is obtained assuming the shape of the 8 B neutrino flux and a total normalization 
10 -2 x < f>( 8 i?) (that is, an overall v e — V e conversion probability P ~ 1%). 


solar magnetic field, this upper bound on the appearance probability implies an upper 
limit on the product of the intrinsic neutrino magnetic moment and the value of the field 
~ fiB < 1CU 21 MeV. For B ~ 10 — 100 kG, we would have n < 10~ n — 10 -12 hb- 

In the opposite case, if spin-flavor precession is indeed at work even at a minor rate, 
the additional flux of antineutrinos could strongly distort the signal spectrum seen at 
KamLAND at energies above 4 MeV and their contribution should be taken into account. 
This is graphically shown in Fig. (6) 
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